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Abstract 

We quantify the ability of i3-Factories to observe supersymmetric contributions to the rare 
decays B — » X s j and B — > X s l + l~ . A global fit to the Wilson coefficients which contribute 
to these decays is performed from Monte Carlo generated data on B(B — » X s ^f) and the 
kinematic distributions associated with the final state lepton pair in B — > X s l + l~ . This fit is 
then compared to supersymmetric predictions. Evaluation of the Wilson coefficients is carried 
out with several different patterns of the superpartner spectrum. We find that i?-Factories will 
be able to probe regions of the SUSY parameter space not accessible to LEP II, the Tevatron, 
and perhaps the LHC. We also employ the recent NLO calculation of the matrix elements for 
B — ► X s j and find the bound m H ± > 300 GeV in two-Higgs-doublet models using present 
data. 



^Work supported by the Department of Energy under contract DE-AC03-76SF00515. 



1 Introduction 



Softly broken supersymmetry (SUSY) is a decoupling theory, thus making it a challenge to search 
for its effects through indirect methods. When competing with standard model tree-level processes, 
the relative shift in an electroweak observable with respect to the standard model (SM) value should 
not be expected to exceed much more than (a 2 /27r) M^/m 2 . In the post-LEP II and Tevatron era, 
if supersymmetry has not been directly observed, then rh > Mjy, and so the relative shift expected 
in standard model observables from virtual supersymmetry is < 0.5%. Although the la bounds on 
sin 2 9w are approaching this level from the analysis of SLC/LEP data, a more statistically significant 
result would be difficult to obtain given the current data sets available. 

Another approach to indirect searches of supersymmetry is to measure observables where su- 
persymmetry and the standard model arise at the same order in perturbation theory. In this case, 
the SUSY contributions do not suffer an extra a/4n reduction compared to the standard model 
amplitudes. The relative ratio between the lowest order standard model amplitudes and supersym- 
metric partner amplitudes could then be 0(1) if m ~ My/- Rare 5-decay measurements provide 
an opportunity for discovering indirect effects of supersymmetry by this second approach. 

Two problems in the past have marred the attempts to use rare B decays as a good probe 
for physics beyond the standard model. The first is limited statistics, or rather the number of 
B mesons available in data sets which can be used to study and obtain good precision on low 
branching fraction modes. The B factories presently under construction, which will collect some 
10 7-8 B mesons per year, will alleviate this issue. The second difficulty is theoretical. Since all the 
processes occur near 5 GeV the uncertainties in the strong interactions can hide even 0(1) effects 
in the electroweak contributions. However, this problem diminishes significantly [|l], ^ [3| with a 
complete program of next-to-leading order (NLO) computations of the QCD corrections to rare B 
decays. For the processes we will consider here, B —* X s 7 and B — > X s £ + £~, these higher order 
calculations have essentially been completed recently 0, [|, |5|. The inclusive decay B — > Xgj has 
been observed by CLEO with a branching fraction of (2.32 ± 0.57 ± 0.35) x 10~ 4 and 95% C.L. 
bounds of 1 x 10~ 4 < B{B — > X s ^) < 4.2 x 10~ 4 . Meanwhile, experiments at e + e~ and hadron 
colliders are closing in on the observation of the exclusive modes B — ► K^£ + £~ with £ = e and fj,, 
respectively. Once this decay is observed, the utilization of the kinematic distributions of the £ + £~ 
pair, such as the lepton pair invariant mass distribution and forward backward asymmetry 0, 0], 



and the tau polarization asymmetry |T0| in B — > X s t + t~ , together with B(B — ► X s ~f) will provide 
a stringent test of the SM. 

The outline of this paper is as follows. In section 2 we calculate the ability of future 5-factories 
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to experimentally determine the magnitude and sign of the relevant Wilson coefficients in the rare 
B decay interaction Hamiltonian using a global fit procedure. We find that the sensitivity for new 
physics will be substantially increased beyond what is currently possible. In section 3 we apply 
these results to supersymmetry and estimate the sensitivity to high supersymmetric mass scales. 
We also reexamine the constraints on the H ± sector from B — > X s 7, in light of the recent NLO 
computations. And in the final section we discuss our conclusions. 



2 Determination of the Wilson coefficients 

The effective field theory for the decays B — > and B — > X s £ + £~, which incorporates the QCD 
corrections, is governed by the Hamiltonian 

n eff 



-r^s^wo*), (i) 

V 2 i=l 



where the Oi are a complete set of renormalized operators of dimension six or less which mediate b — > 
s transitions. These operators are catalogued in, e.g., Ref. [O]. The Cj represent the corresponding 



Wilson coefficients which are evaluated perturbatively at the electroweak scale where the matching 
conditions are imposed and then evolved down to the renormalization scale fi ~ m&. We note that 
the magnetic and chromomagnetic dipole operators, Oy t %, contain explicit mass factors which must 
also be renormalized. 

For B — ► X s [ + [~ this formalism leads to the physical decay amplitude (neglecting the strange 
quark mass) 

m{b - x s i + n = ^^v tb v* \c e 9 ff s Llll b L hn + c I0 ^w 7 v 

(2) 



7T 



-2C e 7 ff m b s L ta^b R !Y£ 



where q 2 represents the momentum transferred to the lepton pair. The expressions for Cj(Mjy) are 
given by the Inami-Lim functions [|l2j. A NLO analysis for this decay has recently been performed^, 



where it is stressed that a scheme independent result can only be obtained by including the leading 
and next-to-leading logarithmic corrections to Cg(fi) while retaining only the leading logarithms in 
the remaining Wilson coefficients. The residual leading \i dependence in Cg(/x) is cancelled by that 
contained in the matrix element of Og. The combination yields an effective value of Cg given by, 

C e 9 ff (s)=CM V (s) + Y(s), (3) 
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with Y(s) being the one- loop matrix element of Og, r](s) represents the single gluon corrections to 
this matrix element, and s = q 2 fm\ is the scaled momentum transferred to the lepton pair. The 
effective value for C? (n) refers to the leading order scheme independent result obtained by Buras 
et al. 0. The corresponding formulae for Cj(/i), Y(s) and r)(s) are collected in Refs. |TT|. The 
operator 1O does not renormalize and hence its corresponding coefficient does not depend on the 
value of fi (except for the \x dependence associated with the definition of the top-quark mass). The 
numerical estimates (in the naive dimensional regularization (NDR) scheme) for these coefficients 
are then (taking m p b ole = 4.87 GeV, m p t hys = 175 GeV, and a s (M z ) = 0.118) 

b 7 if 1 — m b +m b ) — U -<J iZ +0.034 J 

C^ = m b - + 2f) = 4.2lj££, (4) 

and 

C ao (//) = -4.55. (5) 

The reduced scale dependence of the NLO versus the LO corrected coefficients is reflected in the 
deviations AC 9 (/x) < ± 10% and AC? ff {fi) « ±20% as \i is varied in the range mfe/2 < \i < 2m b . 
We find that the coefficients are much less sensitive to the values of the remaining input parameters, 
with AC 9 (m b ), AC% ff {m b ) < 3%, varying a s (M z ) = 0.118±0.003 @|, and mf ys = 175±6 GeV 
P^H corresponding to m t (mt) = 166 ±6 GeV. The resulting inclusive branching fractions (which are 
computed by scaling the width for B — > X s £ + £~ to that for B semi-leptonic decay) are found to be 
(6.25±5;gt) x 10- 6 , (5.73±g;?l) x 10" 6 , and (3.24±g:|J) x 10~ 7 for £ = e,/i, and r, respectively, taking 
into account the above input parameter ranges, as well as B s i = B(B — > Xlv) = (10.23±0.39)% ||16|| . 
and m c /m b = 0.29 ± 0.02 0, [13|]. There are also long distance resonance contributions to B — > 
X s £ + £~, arising from B — > K^'if)"} — > K^£ + £~. These appear as an effective (sl7 m 6l)(£7 m £) 
interaction and are incorporated into Cg^ via the modification Y(s) Y'(s) = V(s) + y res (s), 
where Y res (s) is given in Ref. ||17|| . These pole contributions lead to a significant interference between 
the dispersive part of the resonance and the short distance contributions. However, suitable cuts on 
the lepton pair mass spectrum can cleanly separate the short distance physics from the resonance 
contributions. 

The operator basis for the decay B — > Xgj contains the first eight operators in the effective 
Hamiltonian of Eq. (jH). The leading logarithmic QCD corrections to the decay width have been 
completely resummed, but lead to a sizeable \x dependence of the branching fraction (as demon- 
strated above with the large value of AC? ), and hence it is essential to include the next-to-leading 
order corrections to reduce the theoretical uncertainty. In this case, the calculation of the pertur- 
bative QCD corrections involves several steps, requiring corrections to both the Wilson coefficients 
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and the matrix elements of the operators in Eq. (Jl|) in order to ensure a scheme independent result. 
For the matrix elements, this includes the QCD bremsstrahlung corrections [1^] b — > 57 + g, and 
the NLO virtual corrections which have recently been completed in both the NDR and 't Hooft- 
Veltman schemes Summing these contributions to the matrix elements and expanding them 
around /1 = m b , one arrives at the decay amplitude 

M{b^ S1 ) = - ^^ D( S1 \0 7 (m b )\b) treei (6) 

with 

D = CV e// (/x) + (cf )eff {^ eff log ^ + Cf . (7) 

Here, the quantities "f$ e ^ are the entries of the effective leading order anomalous dimension matrix, 
and the are computed in Greub et al. 0, for i — 2, 7, 8. The first term in Eq. ^, Cy (//), 
must be computed at NLO precision, while it is consistent to use the leading order values of the 
other coefficients. The explicit logarithms a s {m b ) \og(m b / /j,) in the equation are cancelled by the /x 
dependence of C^ e ^(p). This feature significantly reduces the scale dependence of the resulting 
branching fraction. The contribution to the inclusive width including these virtual corrections is 
then 

m b,pole^ 2 F a em \ VtbV t * s | 2 



Tnlo(B - X sl ) = "•^" f 3 J- ttl F\D\> , (8) 

where the factor F = rnl(m b )/ml pole = 1 — 8a s (m b )/3ir arises from the mass factor present in 
the magnetic dipole operator. This should be compared to the familiar leading order result (which 
omits the virtual corrections to (0 7 )) 

r{B - x sl ) = m ^ff Qem |^| 2 |c?^M| 2 . (9) 

For the Wilson coefficients, the NLO result entails the computation of the 0(a s ) terms in 
the matching conditions, and the renormalization group evolution of the Cj(/i) must be computed 
using the 0(a 2 s ) anomalous dimension matrix. The former step has been completed^], but the 
latter step is quite difficult since some entries in the matrix have to be extracted from three-loop 
diagrams. Nonetheless, preliminary NLO results for these anomalous dimensions have recently been 
reported 0, with the conclusion being that in the NDR scheme the NLO correction to C% (/j) is 
small. Therefore, a good approximation for the inclusive width is obtained by employing the 
leading order expression for Cj (/i), with the understanding that this introduces a small inherent 
uncertainty in the calculation. 

The total inclusive width is then given by the sum of the virtual and bremsstrahlung corrections, 
T(B — > Xgj) = Y mrt + T brems , where Y brems is given in Greub et al. |18|, 0], and the branching fraction 
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Figure 1: The branching ratio of B — > X s 7 versus m t . The dashed lines indicate the error in 
the branching ratio if we fix = mj and vary all the other parameters over their allowed ranges: 
a s (M z ) = 0.118 ± 0.003, B sl = 10.23 ± 0.39%, and m c /m b = 0.29 ± 0.02. The solid lines indicate 
the error for m b /2 < fi < 2m b and all other parameters fixed to their central values. 

is calculated by scaling to the semi-leptonic decay rate. The leading order power corrections in the 
heavy quark expansion are identical for B — > X s ~f and B — > X£v, and hence cancel in the ratio |20[] . 
This allows us to approximate T(B — > Xg'j) with the perturbatively calculable free quark decay 
rate. For mf ys = 175 ± 6 GeV, m b /2 < /i < 2m b , a s = 0.118 ± 0.003, B sl = (10.23 ± 0.39)%, and 
m c / m b — 0.29 ± 0.02, we find the branching fraction 

B(B -> X s7 ) = (3.25 ± 0.30 ± 0.40) x 10~ 4 , (10) 

where the first error corresponds to the combined uncertainty associated with the value of mt and 
/i, and the second error represents the uncertainty from the other parameters. This is well within 
the range observed by CLEO. In Fig. |I] we display our results for B(B — > X s , y) as a function of 
the top mass. The dashed lines indicate the error in the branching ratio if we fix \x = m b and 
vary all the other parameters over their allowed ranges given above. The solid lines indicate the 
error for m b /2 < fi < 2m b with all other parameters fixed to their central values. This visually 
demonstrates that the error in the theoretical calculation of B — > X^ is not overwhelmed by the 
scale uncertainty; other uncertainties are now comparable. 

Measurements of B(B — > X^) alone constrain the magnitude, but not the sign, ofC% ff (/i). We 
can write the coefficients at the matching scale in the form Ci(M w ) = Cf M (M w )+C™ ew (M w ), where 
C™ ew (My/) clearly represents the contributions from new interactions. Due to operator mixing, 
B — > Xg'j then limits the possible values for Cf ew {Mw) for % — 7, 8. These bounds are summarized 
in Fig. 0. Here, the solid bands correspond to the constraints obtained from the current CLEO 
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Figure 2: Bounds on the contributions from new physics to C^g. The region allowed by the CLEO 
data corresponds to the area inside the solid diagonal bands. The dashed bands represent the 
constraints when the renormalization scale is set to /i = m^. The diamond at the position (0,0) 
represents the standard model. 



measurement, taking into account the variation of the renormalization scale m&/2 < ft < 2rrib, as 
well as the allowed ranges of the other input parameters. The dashed bands represent the constraints 
when the scale is fixed to ft = m^- We note that large values of C$ ew (Mw) are allowed even in the 
region where Cj ew (Mw) — 0. Experimental bounds on the decay b — ► sg are needed to constraint 

Measurement of the kinematic distributions associated with the final state lepton pair in B — > 
X s £ + £~ as well as the rate for B —>■ X s ^y allows for the determination of the sign and magnitude of all 
the Wilson coefficients for the contributing operators in a model independent fashion || 10|. Here, 
we perform a Monte Carlo analysis in order to ascertain how much quantitative information will 



be obtainable at future 5-factories. We improve upon our previous study[|T0j by implementing the 
NLO computations for these decays and by examining the luminosity dependence of the resulting 
global fits. For the process B — > X s £ + £~, we consider the lepton pair invariant mass distribution and 
forward-backward asymmetry for £ = e, fi, r, and the tau polarization asymmetry for B — > X s t + t~ . 
We note that the asymmetries have the form A{s) ~ C w (ReCg^ fi(s) + Cy" f 2 (s)), and hence are 
sensitive probes of the Wilson coefficients. We generate "data," assuming the SM is realized, by 
dividing the lepton pair invariant mass spectrum into nine bins. Six of the bins are taken to be in 
the low dilepton invariant mass region below the J/ip resonance (in order to take advantage of the 
larger statistics), with 0.02 < s < 0.32 and a bin width of As = 0.05. We have also cut out the 
region near q 2 = in order to remove the photon pole. The high M^+i- region above the ip' pole 
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is divided into three bins, corresponding to 0.6 < s < 0.7, 0.7 < s < 0.8, and 0.8 < s < 1.0. The 
number of events per bin is calculated as 



ds 

and the average value of the asymmetries in each bin is 



£ f s ™«* dT(B -> , , 

{A) bin = — A — — -i '-ds. (12) 



Nbin Jsmin d§ 

We statistically fluctuate the "data" using a normalized Gaussian distributed random number pro- 



cedure, where the statistical errors are given by 5N = y/N and 5 A = y (1 — A 2 )/iV. We expect the 
errors in each bin to be statistics dominated and hence we neglect any possible systematic errors. 
For B — > X s 7, we again statistically fluctuate the "data" for the inclusive rate. However, in this 
case, the statistical precision will eclipse the possible systematic and theoretical accuracy. We thus 
assume a flat 10% error in the determination of the branching fraction in anticipation of future 
theoretical and experimental improvements. A three dimensional x 2 fit to the coefficients CV^io (/•*) 
is performed, employing the usual prescription 




X?=E ~ TFi , (13) 



for each observable quantity Qi. We repeat this procedure for three values of the integrated luminos- 
ity, 3 x 10 7 , 10 8 , and 5 x 10 8 BB pairs, corresponding to the expected e + e~ 5-factory luminosities 
of one year at design, one year at an upgraded accelerator, and the total accumulated luminosity 
at the end of the programs. Hadron colliders will, of course, also contribute to this program, but it 
is more difficult to assess their potential systematic and statistical weights without further study. 

The 95% C.L. allowed regions as projected onto the Cg(fi) — Cio(/i) and (/i) — Cio(n) planes 
are depicted in Figs. |3|(a-b), where the diamond represents the expectations in the SM. We see that 
the determinations are relatively poor for 3 x 10 7 BB pairs and that higher statistics are required 
in order to focus on regions centered around the SM. Clearly, C 9 and C w are highly correlated, 
whereas C7 and C w are not. We see that the sign, as well as the magnitude, of all the coefficients 
including can now be determined. 

For the remainder of this paper, we analyze the supersymmetric contributions to the Wilson 
coefficients in terms of the quantities 

_ cr\M W ) _ cr w (M W ) 
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Figure 3: The 95% C.L. projections in the (a) Cg — Cio and (b) C7 — C10 planes, where the allowed 
regions lie inside of the contours. The solid, dashed, and dotted contours correspond to 3 x 10 7 , 
10 8 , and 5 x 10 8 BB pairs. The SM prediction is labeled by the diamond. 



where C- usy (Mw) includes the full standard model plus superpartner contributions. Ri is meant to 
indicate a relative fraction difference from the standard model value. It is most convenient to define 
these ratios at the W scale to avoid the added complication of the renormalization group evolution 
to the low scale. 



3 Expectations in supersymmetry 

Supersymmetry has many potential sources for flavor violation. The flavor mixing angles among 
the squarks are a priori separate from the CKM angles of the standard model quarks. If we 
allow the super-CKM angles to be arbitrary then we find, for example, that the relative SUSY 
versus SM amplitudes for b — > s are (Mw/m) n (V t *Vtb/V t *V t b)- \V t *V t b\ — 0.04, and so allowing the 
V angles to be oriented randomly with respect to the CKM angles constitutes a flavor problem 
for supersymmetry if m is near the weak scale. Natural solutions to this problem exist. One 
solution is alignment [^1] of the super-CKM matrices with the quark matrices. Another natural 
solution is universality induced by gauge mediated supersymmetry breaking p2| . In the minimal 
model of gauge mediated supersymmetry breaking [Q , the b — > 57 decay is currently not a strong 
constraint on the spectrum, but will show deviations from the standard model at the 5-factory if 



m-± < 350 GeV. We adopt the viewpoint in this paper that flavor-blind (diagonal) soft terms [24 
at the high scale are the phenomenological source for the soft scalar masses at the high scale, and 
that the CKM angles are the only relevant flavor violating sources. We build on other studies of 
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super symmetry effects on rare B decays |25|, ||, |26 

The spectroscopy of the supersymmetric states is model dependent. We will analyze two possi- 
bilities. The first possibility is that all the supersymmetric states follow from common scalar mass 
at high scale and common gaugino mass at the high scale. This is the familiar minimal supergravity 
model. The second possibility is to relax the condition of common scalar masses at the high scale 
and allow them to take on more uncorrelated values at the low scale while still preserving gauge 
invariance. 

We begin by searching over the full parameter space of minimal supergravity model. We use 
the words "minimal supergravity" as an idiom to indicate that we generate these models by 
applying common soft scalar masses and common gaugino masses at the boundary scale. The tri- 
scalar A terms are also an input at the high scale and are universal. The radiative electroweak 
symmetry breaking conditions yield the B and /z 2 terms as output, with a sign(/x) ambiguity left 
over as an input parameter. (Here \i refers to the Higgsino mixing parameter.) We also choose tan j3 
and restrict it to a range which will yield perturbative Yukawa couplings up to the GUT scale. 

We have generated thousands of solution according to the above procedure. The ranges of our 
input parameters are < m < 500 GeV, 50 < m 1( / 2 < 250 GeV, — 3 < A /m < 3, 2 < tan/5 < 50, 
and we have taken m^ hys = 175 GeV. Each supersymmetric solution is kept only if it is not in 
violation with present constraints from SLC/LEP and Tevatron direct sparticle production limits. 
For each of these remaining solutions we now calculate -R7-10 F%fl . Expansions of Ci(M w ) with the 
assumption of approximate universality are given in the appendix. 

First, we present a scatter plot of Rj vs. R$ in Fig. [|; we remind the reader that these quantities 
are evaluated at the electroweak scale. Again, each point in the scatter plot is derived from the 
minimal supergravity model with different initial conditions. Also, each point is consistent with all 
collider bounds and is out of reach of LEPII. The first thing to note from the figure is that large 
values of R 7 and R$ are generated, and the R 7 and R$ values are very strongly correlated. The 
diagonal bands represent the bounds on the Wilson coefficients from the observation of B — > X s ^ 
as determined in the previous section. We note that these bands appear to be straight here as they 
correspond to a small region of Fig. |^. We see that the current CLEO data already places significant 
restrictions on the supersymmetric parameter space. Further constraints will be obtainable once 
a 10% measurement of B(B — > A s 7) is made, and the sign of C-j is determined from a global fit 
described in the previous section. In this case, if no deviations from the SM are observed, the 
supersymmetric contributions will be restricted to lie in the dashed band. It is clear that these 
processes can explore vast regions of the supersymmetric parameter space. In fact, it is possible 
that spectacularly large deviations in rare B decays could be manifest at B factories, while collider 
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Figure 4: Parameter space scatter plot of Rj vs. Rs in minimal supergravity model. The allowed 
region from CLEO data, as obtained in Fig. 2, lies inside the 2 sets of solid diagonal bands. The 
dashed band represents the potential 10% measurement from the previously described global fit to 
the coefficients. 

experiments would not detect a hint of new physics. 

The large effects in Rj and Rs are coming from models with |//| < 400 GeV as can be seen in 
Fig. This is because light charged Higgsinos, or rather light charginos with a large Higgsino 
fraction, are required in order to yield a large effect on the Wilson coefficients. Later in this section 
we will demonstrate this requirement more carefully by expanding the supersymmetric contributions 
in the Higgsino limit. 

In Fig. |6| the correlation between Rg and Rio is plotted using the same supersymmetric parameter 
space. We see that R 9 is always positive since charged Higgs and chargino contributions always add 
constructively. We see that the values of Rg and R w are bounded by about 0.04, a small number 
compared to the range for R7. The main reason for these smaller values is the dependence on the 
bottom Yukawa A& oc 1/ cos/3. R-j has a contribution directly dependent on this 1/ cos/5 Yukawa 
enhancement, and the other multiplicative terms associated with this Yukawa are the large top 
Yukawa and a large kinematic loop factor. R 9 and Rio do not have such factors due the chirality 
structure of these operators and the requirement that leptons and sleptons only couple via SU(2) 
and U(l) gauge couplings. These factors, along with the correlations between the mass spectra 
dictated by minimal supergravity relations, render the minimal supergravity contributions to -R9.10 
essentially unobservable. The solid lines in this figure correspond to the 95% C.L. bounds obtainable 
with high integrated luminosity (5 x 10 8 BB pairs) at B factories from the global fit performed 
in the previous section. If large deviations in -Rg^o are observed, then, of course, that would be 
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Figure 5: Parameter space scatter plot of R 7 vs. ji in minimal supergravity model. 




Figure 6: Parameter space scatter plot of Rg vs. R\o in minimal supergravity model. The global 
fit to the coefficients obtained in Fig. 3 with 5 x 10 8 BB pairs corresponds to the region inside the 
diagonal bands. 



11 




Figure 7: The (a) differential branching fraction, (b) lepton pair forward-backward asymmetry, 
and (c) tau polarization asymmetry as a function of the scaled momentum transfer to the lepton 
pair, s = q 2 /ml. The dotted curves represent the SM prediction, while the dashed and solid bands 
correspond to the maximal potential deviations due to supersymmetric contributions for different 
regions of the parameter space as described in the text. In some cases the dashed line overlaps with 
the solid line. 

an indication that the minimal model presently under discussion is not the correct description of 
nature. Later in this section we will find that even when the mass correlations of the minimal 
supergravity model are lifted, the quantities -Rg,io still cannot be large. 

We next examine the effects of the minimal supergravity model on the kinematic distributions 
for B — » X s £ + £~ . Using our generated models as input, we calculate the maximal deviations from 
the SM for the Mp+i- distribution, lepton pair forward-backward asymmetry, and tau polarization 
asymmetry. These are displayed in Figs. 0( a " c )- Here, the dotted line corresponds to the SM 
prediction, while the dashed (solid) bands represent the maximal possible deviations due to points 
in the supersymmetric parameter space which are within (outside) the expected reach of LEPII. 
We have also demanded consistency with the present CLEO data on B — > X s 7. We see from the 
figure that, generally, larger deviations are expected for the asymmetries than for the branching 
fraction, and that constraints from LEPII on supersymmetry greatly affect the magnitude of these 
potential deviations. Since the SUSY contributions to Rg t io are so small these deviations are mainly 
due to ^7,8- We find that much larger effects in these distributions are possible if the constraint 
from radiative B decays is not taken into account. 

We now adopt a more phenomenological approach. The maximal effects for the parameters 
Ri can be estimated for a superparticle spectrum independent of these high scale assumptions. 
However, we still maintain the assumption that CKM angles alone constitute the sole source of 
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Figure 8: The maximum value of Rg = Cg Usy (Mw) / Cg M (Mw) — 1 achievable for general super- 
symmetric models. The top solid line comes from t — H ± contribution and is displayed versus the 
mass. The bottom solid line is from tj — xf contribution with tan (3 = 1 and is shown versus the 
xf mass. The dashed line is the ti — xf contribution with tan/? = 2. The other mass parameters 
which are not plotted are chosen to be just above LEPII and Tevatron's reach. 

flavor violations in the full supersymmetric lagrangian. We will focus on the region tan f3 < 30 since 
enormous effects are possible in the large tan/3 limit; later on we will discuss the large tan/3 limit 
more carefully. The most important features which result in large effects are a light t\ state present 
in the spectrum and at least one light chargino state. For the dipole moment operators a light 
Higgsino is most important. A pure higgsino and/or pure gaugino state have less of an effect than 
two mixed states when searching for maximal effects in Cg and C\$. In fact, we have found that 
M 2 ~ 2/x is optimal. 

Fig. and |^ display the maximum contribution to Rg y ig = Cg^g (Mw) / Cg%(Mw) — 1 versus 
an applicable SUSY mass scale. The other masses which are not plotted (ti, II, etc.) are chosen 
to be just above the reach of LEPII or the Tevatron, whichever gives better bounds. 

The maximum size of i2g,io is much larger than what was allowed in the minimal supergravity 
model. The reason for this is the lifted restriction on mass correlations. Light sleptons, sneutrinos, 
charginos, and stops are allowed simultaneously with mixing angles giving the maximal contribution 
to the .Rj's. However, we find that the maximum allowed values for -Rg^o are still much less than 
unity. Earlier we determined that B factory data would be sensitive to ARg > 0.3 and Ai? 10 > 0.08 
at the highest luminosities, and so the largest SUSY effect would give a 1 — 2cr signal in -R^io, 
hardly enough to be a compelling indication of physics beyond the standard model. If, on the 
other hand, much larger deviations of -Rg,io are found in the data, it could mean the assumption of 
only CKM angles allowed in the supersymmetric lagrangian is inaccurate, or it could indicate that 
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Figure 9: The maximum value of Rio = C^ sy (Mw) I Cf M (Mw) — 1 achievable for general super- 
symmetric models. The top solid line comes from t — H ± contribution and is shown versus the H ± 
mass. The bottom solid line is from U — xf contribution with tan (3=1 and is plotted versus the 
xf mass. The dashed line is the ti — xf contribution with tan/? = 2. The other mass parameters 
which are not presented are chosen to be just above LEPII and Tevatron's reach. 



minimal supersymmetry is not the source of physics beyond the standard model. 

It should be remembered that even though it appears difficult to resolve the SUSY contributions 
to the coefficients Cg_io, the B — > X s l + l~ decay rate and distributions can still demonstrate large 
deviations from the standard model induced largely by the SUSY corrections to C-jiM^/). The 
global fit using all the rare B decay data is still necessary in this circumstance since it will enable a 



precise determination in which band in the Cj vs. Cs plot we reside. Furthermore, some ideas ]29 
of physics beyond the standard model predict small corrections to B — > X s 7 and large deviations 
in B — > X s l + l~, motivating again the procedure of doing a global fit to all the rare B decay data. 

Given the sensitivity of all the observables it is instructive to narrow the focus to Cr(Mw)- In 
the minimal supergravity models, the scalars obtain dependence on the gaugino masses through 
the renormalization group running. However, the gaugino masses do not get substantial scalar 
mass contributions to their masses. This tends to separate the mass scale for the scalars far from 
the gauginos. The separation is especially important between electroweak gauginos and strongly 
interacting squarks. Neglecting the D-term contributions the squarks have masses given roughly by 
w? q « rriQ + 7my 2 and the weak gaugino has mass ~ 0.8mi/2 where m is the common scalar 
mass and m\ii is the common gaugino mass at the high scale. From these equations it is easy to 
see that the squark masses are much heavier than the weak gaugino mass for any given m 1 / 2 arid 
mo. 

When all the squark masses are very heavy (much heavier than the charginos, for example) then 
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the SUSY contributions to CV(Mw) decouple. However, one eigenvalue of the stop squark mass 
matrix might be much lighter than the other squarks. The large top Yukawa tends to drive much 
lower than the other <1l,r- Also, the stop squark has a mixing term proportional to m t (A t — \i cot /3). 
Since this mixing is proportional to the top mass, it can be substantial. Mixings in any positive 
definite matrix will push the lightest eigenvalue lower and the heaviest eigenvalue higher. These 
two effects tend to push the lightest stop eigenvalue well below the other squarks. In fact, a highly 
mixed, light stop squark is generic in these theories. For a large supersymmetric contribution to 



Rb |30] and/or the ability to achieve successful baryogenesis in the early universe |pl |, a light in 
is needed. We would therefore like to present results on Cj(Mw) in the limit of one light squark, 
namely the ti, and light charginos. We allow the t\ to have arbitrary components of it and tn since 
cross terms can become very important. This is especially noteworthy in the high tan/3 limit as we 
will discuss below. 

In this limit of one light top squark, we can expand the chargino-stop contribution to C 7 {Mw) 

as 



SC 7 (M W ) = W^f {^fn^) +2X^J(x l )\ , (15) 

• 771 ± i mi) ) 

1 Xi 



where 



Af = -TnVa+T, 



V i2 m t 



12" 



X j 
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Ui2 rn h 1 
y/2 M w cos p ' 



2 

m~ 

2 

m ± 



5-7x 2x- 3x 2 , 



_ , N 7 - 5a; - 8x 2 2x - 3x 2 , 
6(1 — x) A (1 — xp 

The matrix is the stop mixing matrix which rotates (IlIr) into {tit 2 ). The matrices JJ^ and 



Vij are the usual chargino mixing matrices |3^j. For the reader's convenience we have plotted the 
functions F(x) and J(x) in Fig. [II]. 

The total contribution to 8C 7 {Mw) above will depend on several combinations of mixing angles: 
T-Q, TnT 12 , Vii?7i2, etc, and cancellations can occur for different signs of /x p4| . Therefore, it is 



instructive to summarize the relative signs of each angle combination in the evaluation of Eq. [H 

sgn(£/i 2 Vn) = -sgn(ii) , 
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Figure 10: The kinematic loop functions F(x) and J(x) necessary to calculate the standard 
model and supersymmetric contributions to C 7 (M W ). J(x) and F(x) asymptote to 5/2 and 7/6 
respectively as x — > 0. 



sgn(L r i 2 V r i 2 ) = +sgn(/i) , 

+, if M 2 tan/? + ii< 0, 
sgn(VnVi2) = -sgn(M 2 tan/3 + /i) , 
sgn(TnTi 2 ) = +, if t R < t L and A t - /x cot /3 < , 
— , if £r < tx and A t — \i cot ft > , 
— , if £x < t/? and A t —\i cot /3 < , 
+, if tx < t R and A 4 — /i cot ft > . 

We are using the Haber-Kane convention for the sign of /i which requires, for example, that +/x be 
in the chargino mass matrix. 

The first case we discuss is the limit where the lightest chargino is a pure Higgsino and the 
lightest stop is pure right-handed: Xi ~ H ± , t% ~ t R . Then Eq. |15| can be written as 



5C 7 {M ] 



w , 



12 m 2 



2 



ti 

2 

m + 

Xi 



(16) 



This limit is roughly applicable in the case where i4 is affected substantially by supersymmetric 



corrections |30], |33|. After LEP II completes its run the charginos will have been probed up to 
about M\y, and so the maximum effect on 5Cj(M\y) that would be possible in this limit (if LEP II 
does not find charginos) corresponds to setting m x ± ~ Mjy. We do this and show the result as a 
function of the t R mass in Fig. [Il| (dashed line). The contribution to Cj(Mw) in this limit is always 



16 




50 100 200 500 1000 2000 

Mass (GeV) 

Figure 11: Contributions to R-j in different limits. The top solid line is the charged H ± /t contri- 
bution versus m#±. The bottom solid line is the xt/ti contribution versus m^± where both the 
chargino and stop are maximally mixed states (equal magnitude mixtures for all states). This line 
was made for /x < 0. The dashed line is the H /£# contribution, and the dotted line is the W /t\ 
contribution (ti oc tn + tz, is a maximally mixed stop mass eigenstate). These two lines are both 
plotted against the xf mass. All lines are for tan/5 = 2 and m t = 175 GeV. We have set all other 
masses to be just above the reach of LEPII. 

positive. Since C 7 M (M W ) is a negative quantity in our convention, then R 7 = 5C 7 (M W ) / C 7 M (M w ) 
is necessarily negative as well. 

In the case where the only light chargino is a pure Wino we find, 

SC 7 (M W ) = T ±<L F . (17) 

6 m w \ m wJ 

If li is very large and tan/3 is not too large, then this contribution will be the largest to SC 7 (Mw)- 
The effects of a light pure Wino are generally small since the Wino couples like #2 rather than the top 
Yukawa, and since generally supersymmetric models do not yield light necessary to couple with 
the Wino. The loop integral F(x) which characterizes the pure gaugino contribution is also small. 
Therefore, these contributions are both coupling suppressed and "loop integral" suppressed. The 
contribution of a pure Wino to R 7 is shown in Fig. [Tl] (dotted line). As expected, this contribution 
is rather small and always negative. 

As we pointed out above, in minimal supergravity models, what we mostly expect is a highly 
mixed t\ state such that all entries in the are filled. By just looking at the limit of pure 
tji some very interesting effects can be missed. For example, an attractive high scale theory is 
supersymmetric 5*0(10). These have been shown to successfully recover many important features 
of the standard model: quark to lepton ratios, CKM angles, etc. One of the generic predictions in 
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these models is that tan (3 must be rather high (~ m t jm^) in order to get the b—r—t mass unification. 
Therefore, it is interesting to focus on contributions to Ci(Mw) which are especially dependent on 
high tan/3. Since Af oc 1/ cos/5 ~ tan/3, we isolate this piece. For simplicity we expand about 
xf ~ H ± , although it is clear that substantial contributions exist even if the chargino is not a pure 
Higgsino. We then obtain 

,2 



5C 7 ~ sgn(//)-— ^T n T 12 tan /? J . (18) 



TnTi 2 tan/3J -^j^- 

Here it is crucial that there be substantial £r and contributions to £i. As argued above, this is 
generic in minimal supergravity theories. This expansion demonstrates that large tan/3 solutions 
(tan/3 > 40) can yield greater than 0(1) contributions to Rj for mass scales of ITeV. Even rather 
low values of tan/3 exhibit enhancements with a light Higgsino and light mixed stop. This is 
demonstrated for tan/3 = 2 in Fig. |ll] (solid line). Furthermore, large contributions are possible 
in both the negative and positive directions of R-j depending on the sign of fi. For example, with 
m h = 250 GeV, \T U T 12 \ = 1/10, m-± = 250 GeV and tan/? = 50 we find that \R 7 (M W )\ > 3. 
Again, we are in a region of parameter space which is not tuned just to give this large effect in 
B — > A s 7, rather we are in a region which is highly motivated by SO( 10) grand unified theories. 



The finite corrections to the b quark mass [35, [36] constitute approximately a 20% correction to the 
b Yukawa coupling when tan {3 > 40. This b Yukawa coupling is implicitly present in \f . Depending 
on the sign of /i this correction can be positive or negative. We don't include these finite b mass 
corrections in our analysis since it requires a detailed knowledge of the sparticle spectrum which we 
are not specifying. 

The chirality structure of the Oi operators allow a large tan /3 enhancement only for the C? t s {My/) 
coefficients. In O7 the bn quark must undergo an helicity flip as long as we neglect m s dependent 
effects. Therefore, all contributions to C7 in the standard model and MSSM must be proportional 
to the 6-quark mass. However, some diagrams with tk/xf loops allow proportionality to the bottom 
Yukawa alone, which yields a m&/ cos/3 enhancement for large tan/3 [|37|]. The mixings between the 
charged Higgs weak eigenstate and goldstone forbid vn^j cos /3 enhancements in the physical charged 
Higgs graphs. Furthermore, the helicity structure of the four-fermion operators O910 forbid large 
tan/3 enhancements. 

We conclude our analysis by examining the charged Higgs contributions to B — > A s 7 alone. It 
is well-known [ 38], |39|j that a H ± boson can contribute significantly to CV,8, but has a smaller effect 
on the coefficients Cg^o; this is also illustrated above in Figs. 8,9,11. The form of the coefficients of 
the magnetic dipole operators in this case are presented in the appendix. From these equations, we 
see that not only do large enhancements occur for small values of tan j3, but more importantly, the 
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Figure 12: Constraints in the charged Higgs mass - tan/? plane from the CLEO bound on B(B 
X s j). The excluded region is that to the left and below the curves. The top line is for m 



phys 



181 GeV and the bottom line is for 



0.52 GeV" 



169 GeV. We also display the restriction tan (3/m H ± > 
which arises from measurements of B 



Xtv as discussed in [40 



coefficients are always larger than those of the SM, independent of the value of tan/3. This leads 
to the familiar bound m#± > 260 GeV obtained from the measurement of BiB — > X^) by CLEO. 
However, this constraint does not make use of the recent NLO calculation of the matrix elements for 
this decay which are discussed in previous sections. We remind the reader that a full NLO calculation 
would also require the higher order matching conditions for the SM and H ± contributions as well 
as the NLO anomalous dimensions for Cy (/i). Nevertheless, we recall that preliminary results on 
the NLO corrections to indicate they are small 0, and a good approximation is obtained 

by employing the leading order expression for Cj^(fi) with the NLO matrix elements. Since this 
drastically reduces the /i dependence of the branching fraction, we would expect the H ± constraints 
to improve. Indeed, we find that the CLEO bound excludes the region to the left and beneath the 
curves in Fig. [T2|. For m1 hys = 169 GeV we see that w,h± > 300 GeV. This is calculated by using the 
same procedure that produced the previous charged Higgs mass bound by CLEO, i.e., all the input 
parameters (e.g., a s , fi, m c /mb, and B(B — > X£u)) are varied over their allowed ranges in order to 
ascertain the most conservative limit. This bound holds in the general two-Higgs-doub let-model II, 
and in supersymmetry if the superpartners are all significantly massive. 
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4 Conclusion 



In this paper we have studied the effects of super symmetry to the FCNC observables concerning 
b —>■ s transitions, and we have seen that deviations from the standard model could be detected 
with supersymmetric masses even at the TeV scale. This is especially true if tan f3 is very high. 
The large tan/3 enhancements in the b — > s processes are unique. Deviations in B — B mixing 
are not as pronounced since the m&/ cos (3 enhancements, which are possible in this case, must be 
compared with m t , whereas in the magnetic dipole b — > s transitions the m&/ cos/3 enhancement 
is to be compared with m&. Furthermore, the uncertainties in the decay constant and bag factor 
make B — B mixing observables slightly less appealing when trying to probe deviations from the 
standard model. 

Therefore, it is possible that the first distinct signs of supersymmetry could come from deviations 
in b — > s transitions. One would like direct confirmation of a possible deviation at the S-factory, 
and collider programs such as an upgraded luminosity Tevatron [M]], and the LHC could provide it. 
At the Tevatron, the cleanest signal for supersymmetry is the 3/ signal coming from leptonic decays 
of pp — > xfxt- ^ the light charginos and neutralinos are mostly gaugino-like then the branching 
fraction into 3/ can be quite high. This is true when the sleptons are lighter than the squarks and 
near Mw, as the slepton mediated decays of the gauginos enhance the final state leptons branching 
fraction. If the charginos and neutralinos contained a substantial Higgsino fraction then the slepton 
mediated exchanges are suppressed by lepton Yukawas to a negligible level, and all the decays must 
proceed through W and Z bosons, and recall that the branching fraction of WZ — > l + l~V < 2% 
(not counting r's.) On the other hand, the b — > s transitions are largest when there is a substantial 
Higgsino fraction in the lightest chargino state. The good news is that the trilepton signal and 
b — ► s decays are somewhat complementary in supersymmetric parameter space. The bad news is 
that confirmation between the two experiments might be difficult. Confirmation could be possible 
later at the LHC through total missing energy signatures, or at the NLC through direct mass and 
mixing angle measurements which would allow SUSY predicted b — > s rates to be compared with 
the data. 

Much of the signal phenomenology at high energy colliders such as the Tevatron and LHC rely 
on missing energy signatures. However, these are only possible with a stable LSP, which is the 
result of a postulated R-parity among the fields. R-parity conservation certainly need not hold in 
the correct theory of nature, and even tiny R-parity violating couplings destroy the missing energy 
signatures. More complicated analyses then must be performed with greatly reduced sensitivity 
to sparticle production. The Tevatron upgrade then would have difficulty exceeding equivalent 
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LEPII bounds on sparticle masses jEJ . The LHC would have a significant search capability beyond 
LEPII, although much reduced compared to search capabilities with R-parity conservation if the 
LSP decays hadronically It appears to be difficult to extract SUSY signals at the LHC for 
superpartners above 1 TeV in this case. Searches for virtual sparticles, such as those we discuss in 
this paper, do not suffer in the presence of R-parity violation, and we have already noted that SUSY 
contributions are resolvable with masses above 1 TeV. In fact, the signal may be enhanced over 
the rate for gauge interactions alone if the R-parity violating couplings are sufficiently large. It is 
therefore possible that the 5-factory could be sensitive to some parts in supersymmetric parameter 
space not accessible at the LHC. 
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Appendix 



In this appendix all the standard model and supersymmetric contributions to the matching condi- 
tions C 7 ^i (M w ) are listed in the limit that all squarks are degenerate except the stop squarks. It 
is also assumed that no sources of flavor violation are allowed other than the CKM angles. 

The dipole moment operators Cr,s(Mw) are already provided in the literature in several places 
for this limit. For completeness we write them down here using the formulas of [44j which are 
normalized according to our definition above: 
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It is convenient to write Cg^M^) as 
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The values for Y and Z are contained in |26| for arbitrary flavor structure and masses. For all 
squarks degenerate except the stop squarks and with only CKM flavor violation then 
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The functional form of each of these contributions is 
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In the above equations 777.3 is the common squark mass. The functions are defined as 
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These L-f unctions are expressed in terms of functions contained in [pH] and are explicitly given 
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